A constrained random-force model for weakly bending semiflexible polymers 



O 



o 



1 

o 
o 



> 

00 

m 
en 

o 



X 




Panayotis Benetatos and Eugene M. Terentjev 
Cavendish Laboratory, University of Cambridge, J. J. Thomson Avenue, Cambridge, CB3 OHE, United Kingdom 

(Dated: July 19, 2011) 

The random- force (Larkin) model of a directed elastic string subject to quenched random forces 
in the transverse directions has been a paradigm in the statistical physics of disordered systems. 
In this brief note, we investigate a modified version of the above model where the total transverse 
force along the polymer contour and the related total torque, in each realization of disorder, vanish. 
We discuss the merits of adding these constraints and show that they leave the qualitative behavior 
in the strong stretching regime unchanged, but they reduce the effects of the random force by 
significant numerical prefactors. We also show that a transverse random force effectively makes the 
filament softer to compression by inducing undulations. We calculate the related linear compression 
coefRcient in both the usual and the constrained random force model. 

PACS numbers: 05.20.-y, 36.20.Ey, 87.15.ad 



The random-force model was introduced by Larkin as 
a short-distance approximation to study the effect of 
a quenched random potential on the Abrikosov lattice 
[l|, |2| . Because of its simplicity (the relevant functional 
integrals are Gaussian), it admits exact solutions and has 
become a paradigm in the physics of disordered elastic 
manifolds '3| . In a recent publication [4] , we used it to an- 
alytically investigate the effect of a quenched disordered 
environment on a strongly stretched wormlike chain. Of 
course, the randomness that biopolymers are subject to 
in the cellular environment is much more complicated 
[5|, ISl , but this model yields some very simple analytical 
results for the weakly bending case. 

The random-force model is mathematically well- 
defined and its analysis is valid. However, its relevance 
to experimental measurements is obscured by the fact 
that the distribution of quenched disorder includes real- 
izations of the random force which have non-zero total 
force along the polymer contour. The total transverse 
force in any given experiment is either zero or non-zero 
and, irrespective of the size of the polymer, these two 
distinct possibilities persist and do not average. We can 
show, using the central limit theorem, that the variance 
of the total transverse force, in the limit of long contour 
length L, scales as ~ L. This is a manifestation of lack 
of self-averaging, and the disorder-averaged value of an 
observable is not expected to coincide with the outcome 
of a single measurement in a long filament. Moreover, we 
should consider two distinct scenarios in a stretching ex- 
periment. In the first, the end points of the polymer are 
fixed in space and absorb (balance) both the total force 
and torque that the polymer feels. A net transverse ran- 
dom force will have a macroscopically manifest effect on 
the conformation of the polymer, as shown in Fig. [TJ It 
will also be experienced by the end-point clamps. The 
same is true for a net torque. In the second scenario, 
the polymer ends are free to equilibrate in the random 
potential and the net force and torque identically van- 
ish. Averaging over all possible random force realizations 



FIG. 1: An unconstrained random linear force density in- 
cludes realizations with a finite net transverse force balanced 
by the clamps. Such a net force yields an overall bend in 
addition to the microscopic undulations. 



could be misleading in the case of vanishing net force or 
torque. These observations motivate us to study a mod- 
ified version of the random-force model which includes 
the constraint of vanishing total force along the polymer 
contour in every single realization of the quenched disor- 
der [Tj. We also consider the additional constraint pre- 
scribed by the condition of mechanical equilibrium that 
the total torque on the polymer in every single realiza- 
tion of the quenched disorder vanishes. The only effect 
of such a random force is the deformation of the poly- 
mer through random undulations. A quenched trans- 
verse random force could also result from the interaction 
of a stretched semiflexible polyampholyte with a trans- 
verse electric field (assuming screened intrapolymer in- 
teractions). It is known that the constraint of global 
neutrality may affect the behavior of a polyampholyte 
[8[. This provides further motivation for studying the 
constrained random force model. 

We condider a weakly bending wormlike chain of total 
contour length L in 1 -I- 1 dimensions parametrized by 
the displacement y{s) perpendicular to the aligning di- 
rection {x), where s is the arc- length along the polymer 
contour. An axial force pair / is applied at its end-points 
in the x direction and a spatially varying random linear 
force density g{s) is applied along the polymer contour 



in the transverse direction. The elastic energy functional 
is given by 
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where k is the bending rigidity related to the persis- 
tence length Lp via k — ^LpkuT. / > corresponds to 
stretching and / < corresponds to compression along 
the filament axis. The first term in the r.h.s. of Eq. ([1]) 
penalizes bending, the second term expresses the interac- 
tion with the transverse random force, and the remain- 
ing two express the interaction with the stretching force 
{—f[x{L) — x{0)]). g{s) acts as quenched disorder. Its 
probability distribution is given by 
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This is Gaussian at any arc-length position s with vari- 
ance Ag, modified by the constraint of zero total force 
which is expressed by the integral over A using the Fourier 
representation of the Dirac ^-function, and the constraint 
of zero total torque which is expressed by the integral over 
H- For the sake of simplicity, we consider hinged-hinged 
boundary conditions: 

y{0) - y{L) = 

dMB)U=o = d^Ms)U=L=0. (3) 

The force-extension relationship is obtained from 
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where (...) and (...) denote thermal and disorder averages, 
respectively. The boundary conditions of Eq. ([3|) allow 
us to to expand y{s) in a series of sines: 



It: 
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We calculate disorder-averaged correlators using the 
replica method as in Ref. [J|: 
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lab being an n x n matrix with all of its elements equal 
to 1, and 
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The first term in the r.h.s. of Eq. ([9]) is related to the 
usual random-force model whereas the remaining terms 
express the contribution of the constraints. 

In the absence of the zero-torque constraint (that is, 
without the second factor in the rhs of Eq. ^), Eqs. ([7]) 
and ([9]) become 
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respectively. 

In the strong stretching regime which is defined by 
/ :^ n/Lp for filaments with L >• Lp or / >• k/L^ for 
filaments with Lp '^ L, we obtain: 
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The last term in the r.h.s. of the above equation is the 
decrease in the end-to-end distance due to the undula- 
tions caused by the transverse random force. In the 
usual (unconstrained) random-force model, the numer- 
ical prefactor of that term is 1/12 [9j. The zero-total- 
force constraint merely reduces this prefactor to 1/24, 
and adding the zero-total-torque constraint further re- 
duces it to 1/30. 

We also investigate the effect of the constraints on the 
puUing-force response of the width of the transverse un- 
dulations (mean square transverse displacement at the 
mid-point) of the stretched wormlike chain. In the strong 
stretching regime, we obtain: 
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In a similar fashion as with the force-extension response, 
the zero-total- force constraint, reduces the effect of the 
random force, which is expressed by the last term, by 



(9) 



a factor of 1/2 (it changes the numerical prefactor from 
1/48 to 1/96). The zero-total-torque constraint does not 
affect the width of transverse fluctuations. We point out, 
however, that these simple results (Eqs. (|12l) and (IT^ ) 
only hold in the strong stretching regime. For arbitrary 
stretching force /, the difference between the ususal and 
the constrained random force model is more complicated 
and depends on the system size (L) as can be seen from 
Eq. (HI). 

A wormlike chain at zero temperature does not yield to 
a small axial c omp ressional force below the critical buck- 
ling value 10|, lul- At any finite temperature, however. 



thermally induced undulations smooth out the buckling 
and give rise to a linear response for small compressional 



forces [12]. A quenched disordered environment induces 



similar undulations and therefore a similar contribution 
to the linear response coefficient. We consider the con- 
strained random-force as a simple model for the quenched 
disordered environment. Assuming / < (compression) 
and Lp ^ L, we obtain the average projected length of 
the filament in the direction of the compressing force to 
leading (linear) order in |/|: 
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The first line Eq. (ITil) expresses the contribution of the 
thermally induced undulations whereas the second re- 
fiects the effect of the transverse random force. Both 
contributions modify the classical Euler behavior of an 
elastic rod before buckling. This second line becomes 
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in the usual (unconstrained) random force model and 
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if we add the zero-total- force constraint. In contrast to 
the non-linear stretching response (Eq. ([12])), where the 
contribution of the random-force undulations decreases 
relative to that of their thermal counterpart as the pulling 
force increases, in the linear compressional response the 
the two contributions enter on equal footing. We point 
out that the thermal and the quenched-disorder contri- 
butions differ qualitatively as far as their dependence on 
the filament parameters (bending stiffness n, and total 
contour length L) is concerned. This difference could in 
principle be useful to probe the existence of non-thermal 
lateral random forces on compressed microtubules in the 
cellular environment. For f — and T = 0, a filament of 
fixed bending stiffness tends to crumple under the load of 
the random forces as its size L increases. This is typical of 
the Larkin model and it remains unaffected by the zero- 
total-force and zero-total-torque constraints. Of course. 



within the context of our weakly-bending approximation 
we can only identify a tendency towards crumpling and 
not an actual transition. 

In this brief note, we have investigated a modifica- 
tion of the random-force model which excludes realiza- 
tions of disorder with a non- vanishing total force and a 
non-vanishing total torque. These are non-self-averaging 
contributions. We have calculated the force-extension 
response of a wormlike chain and the width of its trans- 
verse undulations in the strong stretching regime using 
this modified model. We have also calculated the linear 
compressional response using both the usual and the con- 
strained random force model. We have shown that the 
constraints leave the calculated behavior qualitatively 
unchanged but they reduce the effect of the random forces 
by significant numerical prefactors. The dependence of 
the disorder-induced contribution on the size of the sys- 
tem (L), which grows in the thermodynamic limit, sug- 
gests that the random force model, even in its modified 
version, is inherently non-self-averaging. It would be in- 
teresting to calculate how the constraints discussed here 
would modify the free energy distribution functions ana- 
lyzed in Ref. |3[. Although in this work we deal with a 
weakly bending wormlike chain, our analysis of the van- 
ishing total force and torque constraints also holds in 
other applications of the random- force model (e.g., flux 
lines in type II superconductors or liquid crystals in ran- 
dom media [13|). 
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